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ABSTRACT 

The aim of this work is to show, in the framework of a simple finite-width disk model, 
that the amount of mass seen through gravitational microlensing measurements in the 
region < R < R Q is consistent with the dynamical mass ascertained from Galaxy 
rotation after subtracting gas contribution. Since microlensing only detects compact 
objects, this result suggests that a non-baryonic mass component may be negligible 
in this region. 

The definitive version is available at 

|http : //dx . doi . org/10 . 1051/0004-6361/201219926 
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1 INTRODUCTION 

Spiral galaxies are customarily regarded as disk-like objects 
imbedded in a spheroidal halo of non-baryonic dark matter 
(NDM). In this picture the rotation of galaxies is governed 
mainly by the gravitational field of NDM halo, at least at 
larger radii. 

But this scenario might not be typical, because ro- 
tation of some of spiral galaxies can be explained with 
a negligible (if any) amount of NDM. For these galaxies, 
a thin disk model approximatio n has p roven to be quite 
satisfactory (|jaIocha et al.1 |200S| . l2010al k There was also 
a successful attempt at explaining high transverse gradi- 
ents i n the Galaxy rotation in the same disk model frame- 
work (Jal ocha et al.ll2~010bl k whereas including a significant 
spheroidal mass component proved to reduce the predicted 
gradient to values that are in consistent with measure ments. 
There is also a recent result l|Moni Bidin et al]|2012l ) show- 
ing that the kinematics of the thick disk stars in the so- 
lar neighborhood is consistent with the visible mass alone, 
and there is no place for a spheroidal distribution of NDM, 
althoug h the assumptions of th is result are the subject of 
debate (iBovv fc Tremaindliorl '). 

The above results suggest that the overall mass distribu- 
tion in Galaxy might be flattened rather than spheroidal, at 
least in the Galaxy interior. In this case, the disk model de- 
scription would be more appropriate to describing the true 
mass distribution in this region than mass models with a 
dominating spheroidal mass component. Further indepen- 
dent arguments that would help for deciding about the prop- 
erties of mass distribution in Galaxy are necessary. To this 



end we apply the gravitational microlensing method as a 
means to determining the mass residing in compact objects, 
and compare the results for radii R < R a with predictions 
of a finite-width disk model. 

It is important that in the same region where mi- 
crolensing results are known, the Galactic rotation curve 
is relatively well determined. Figure [1] shows rotation data 
for the Galaxy from various independ e nt measurements 



( Burton fc Gordon! Il978l; 
Demers fc Battinellill2007l; 



Clemens! 119851; iFich et al.1 119891: 
Blitz et al.lll982l ; lHonma fc Sofud 



Il997[). together with two model rotation cur ve fits obtained 
in l|Sofue et"aill2009l ) and (|Sofue et al.lll999l k which we use 
in our calculations and refer to as A and B, respectively. 
They are used, because they agree well with the data points 
in the region of interest. Rotation curve A was obtained 
assuming the reference velocity = 200km/s. Rotation 
curve B assumes = 220km/s. Our results are based on 
these two curves, although in sec !2.3l in order to estimate 
the influence of the uncertainty in the determination of the 
parameter Oo on our results conce rning the optica l depth, 
we also use a rotation curve from (McMillan 2011J) that is 
based on the terminal rotation curve data and a high value 
Oo = 239km/s preferred by the author. 

The surface mass density in thin disk model can be 
obtained with the help of the following integral mapping 
v(R) - the tangential component of Galaxy rotation in the 
disk plane - to the surface mass density <j(R) in the disk 
plane (representing a column mass density): 
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Figure 1. Rotation of Galaxy. Unified rotation velocity measure- 
ments [points] (see the text for re ferences) and two model rotation 
curves A \thick lin e] dSofue et al.ll2009f) and B [dot-dashed line] 
jSofue et al.lll999t ) used to obtain surface mass density. 
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Here, 7? is the radial variable in the disk plane, K and E 
are complete elliptic integrals of the first and second kindQ 
= f^pf and i is a dimensionless integration variable. 
Mathematical details are given in Appendix [X] In studying 
galaxies with the use of the expression ([TJ , it is assumed that 
a galaxy is a flattened disk-like object rather than spheroidal 
and that the motion of matter is predominantly circular. 
The central bulge can also be described in this model by 
the appropriate substitute column mass density. One should 
recall that the relation between the mass distribution and 
the rotation law is qualitatively different from that valid 
under spherical symmetry. It is nonlocal in the sense that the 
total rotation field is needed to determine the surface mass 
density at a given radius. Cutting off integration introduces 
some error, however it is usually not significant for rad ii 
smaller than 0.6 of the cutoff radius (Bra tek et al.l l2008). 
This condition is satisfied in the region of our interest. The 
resulting surface mass densities corresponding to rotation 
curves A and B are shown in figure [2] 



2 THE MICROLENSING METHOD 

In contrast to methods based on analyzing the Galactic rota- 
tion or the motions of halo objects, the microlensing method 



1 Since there are two conventions for the argument of elliptic 
inte grals met in the literatu re, we give the definitions which we 
use jGradshtci n et alj|2007t) : 
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Figure 2. Surface mass densities in the thin disk model corre- 
sponding to rotation curves A [solid line] and B [dot-dashed line]. 



gives us the opportunity of measuring Galaxy mass inde- 
pendently of Galaxy dynamics. The meth od was develope d 
and described in de tail by several authors (|Paczvnskilll986l ) , 
l|Paczvnskilll996l l. ijSchneider et all 120061 ). To estimate the 
integrated mass distribution along various lines of sight, the 
method employs the phenomenon of gravitational light de- 
flection of distant sources of light in the field of point-like 
masses distributed in space. A directly measured quantity in 
this method is the optical depth. Before going to the main 
analysis, let us discuss this important notion in more detail. 



2.1 Optical depth along a line of sight 

Suppose a point-like lens of mass M is located exactly on 
a line of sight between the source of light and the observer. 
This line plays the role of the optical axis. Then, the image 
of the source forms the Einstein ring of some angular ra- 
dius 6e as perceived by the observer and concentric with 
the optical axis. Let Dl and Ds be the distances along 
the optical axis between the observer and the lens and be- 
tween the observer and the source, respectively. Expressed 
in terms of 0e and the gravitational deflection angle oie, the 
lens equation in the weak field approximation reduces in this 
particular configuration to Qe Ds = o.e{Ds — Dl)- Taking 
into account that ole = with Re = OeDl, we obtain 
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The optical depth is by definition the probability of find- 
ing a compact object (a lens) within its Einstein radius on 
some plane perpendicular to the line of sight. In a given vol- 
ume element dSdDL, large enough to contain many lenses 
of various masses, the corresponding probability is a sum 
ViTiR Ei dDL over all lenses of mass Mi, where Vi is the 
number density of lenses of mass Mi. The ratio R%./Mi is 
independent of Mi, and ^ i ^M, is simply the volume mass 
density of compact objects (lenses) p co . Integration over Dl 
along a given optical axis therefore gives us 
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This is the (integrated) optical depth along the line of sight. 

For practical reasons, the passage of a lens in front of a 
source of light is assumed to have occurred when the mag- 
nification of the source exceeds some threshold value. The 
optical depth is measured by means of counting the number 
of such events for millions of stars located near the Galactic 
center. A compr ehensive descr iption of such a measurement 
can be found in lMonied (|201Ch . 



2.2 Application in the framework of a finite-width 
disk model 

A thin disk model predicts some column density a from 
Galaxy rotation data with the help of the formula discussed 
earlier. To compare the disk model predictions with the 
optical depth measurements, an appropriate volume mass 
density entering the integral defining r is needed. In the 
disk approximation, a can be regarded as a column mass 
density assigned to a volume mass density p. This identi- 
fication is also frequently made in modeling galactic disks. 
The column mass density is defined by means of the inte- 
gral a(R) = p(R, Z) dZ, where R, Z are Galacto-centric 
cylindrical coordinates and the axial symmetry has been as- 
sumed. 

The integration in r is carried out along a line of 
sight joining the observer located at Xq = [7? o ,0,0] 
and a source located at Xq, where X® = Xq + (1 + 
x)Ro [— cos &cos I, — cos 6sinZ,sin &]. Here, 6 is the Galactic 
latitude, / the Galactic longitude, and (1 + x)-R° is the dis- 
tance of the source of light from the Sun (\ is a dimension- 
less position parameter). A convenient parametric descrip- 
tion of the line of sight is X(s) = Xq + s(X® — Xq) (then 
Dl{s) = s(l + x)Ro and Ds = (l + x)R°)- On assuming the 
standard exponential profile p(R, Z) = p(R, 0) e~' z ' // ' 1 with 
a fixed scale height h, the line integral r becomes a function 
of 6, I and \- 



that the disk model gives the surface density describing the 
gravitating masses, regardless of their nature, whereas the 
microlensing events only concern compact objects. Accord- 
ingly, in making comparisons, the contribution from con- 
tinuous media such as gas should be subtracted from the 
density. 

For estimating the gas contribution we used the 
neutral and molecular hydro gen density profiles ob- 
tained by iMisiriotis et all (|2006l ) from COBE/DIRBE and 
COBE/FIRAS observations. The column mass density of 
gas is then found to be approximately four third of that 
for hydrogen. It includes the contribution from helium. The 
result of this procedure is shown in figure [3] 
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Figure 3. The ratio of the hydrogen surface mass density ajj 
relative to the total surface density accounting for the Galaxy 
rotation in the thin disk model. A factor 4/3 has been included 
to account for the helium abundance. 



The relation 77,(6) is plotted in figure [4] for several values 
of b. As is seen, the predicted optical depth has a correct 
(2)order of magnitude. Its dependence on the scale height is 
significantly weaker at higher latitudes. 



a (r o \JI + (1 + x) s cos b [(1 + x)s cos 6—2 cos I ]^ s(I — s)ds 

where a is a known function. 

The observed sources of light are located in the vicin- 
ity of the Galactic center. The measurements we use are 
presented as pairs (t, b) for various b with some errors At. 
They are taken from (|Moniejl2010h as averages over a do- 
main of I around Galactic center. All of the fields used 
are placed within the a rea presented in the figure 14 in 
lHamadache et all (|2006T ). although the most important re- 
gion is I £ (—5°, 5°). The precise value of \ i s n °t known 
and its uncertainty must somehow be taken into account 
in the model of r. Either way, in a first approximation one 
can assume 1 = and x = 0, in which case ([2]) reduces to 
77,(0, &, 0) w 77,(6), where 
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77,(6) = 
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s(l-s)ds. (3) 



In 17,(6) we used the approximations cos 6 ~ 1 and sin 6 m 6, 
since |6| < 0.1 radians for the observed sources of light. 

Before using the above integrals, a remark should be 
made concerning the appropriate choice for a. We recall 



h[pc] 

Figure 4. The optical depth r h (b) in disk model as a function of 
the disk scale height h shown for various lines of sight, labeled by 
the corresponding latitudes b (in degrees). 



To find a better model of the optical depth, the uncertainty 
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in x in the integral ([2]) and an averaging procedure of the 
observational data over I must be simulated. To this end a 
Monte Carlo method can be used. 



2.2.1 A Monte Carlo method 

For a given value of b, a collection of lines of sight (called 
ensemble) is chosen randomly with some probability distri- 
bution for various pairs (Z,x)i so that X® is always in some 
region in the vicinity of the Galactic center. Then t (I, b, \) 
is calculated for every element of the ensemble, and next, 
the average value over the ensemble and the respective dis- 
persion is determined. In our simulation we assumed a uni- 
form distribution for \ an d l- \x\ < 0.125 (that is, zblkpc), 
\l\ < 0.017 (that is, ±l°fl We also used an auxiliary uniform 
random variable E £ (0,<r(0)). At a fixed 6 we choose those 
from among the random triples (x, I, E) for which E < a(R). 
As a result, the distribution of X® is weighted by o. 



2.3 The results 

The frequently used value for the scale height is h = 325 pc. 
This standard choice is consistent with the sca le height pa- 
ramet er for M-dwarfs of 320 ±50 pc obtained bv lGould et al.l 
l|l997h . 

The optical depth calculated in accord with equation 
((3| is plotted in figure [5] for h = 325 pc and compared 
with the optical depths measure d by sev eral collaborations : 
MACH O dPopowski et ai]|2005t> . OGLE (ISumi et al.ll200fj) 
ERO S (jHamadache et al.1 12000 ) . and MOA (|Sumi fc at all 
I2003TI . These data were collected and discussed in a review 



l Moniej[2oTol 'l. There is a subclass in the data consisting 
of bright stars. Separation of bright sources allows reducing 
the blending effect, w hich is a seri ou s source of uncer tainty 
in the optical depth (|Alardl fl997t ). (|Smith et al.l l2007f ). It 
is seen that the optical depth predicted in the framework 
of our model fits the data quite well when the bright stars 
are considered. In figure [5] the result of the Monte Carlo 
procedure described earlier was shown and compared with 
the result shown in figure [S] The simulation reproduces the 
optical depths consistently with measurements. The Monte 
Carlo standard deviation also takes some information about 
the expected uncertainty into account in the distribution of 
sources around the Galactic center (as a function of lon- 
gitude and distance to the source Ds). It is worth noting 
that this uncertainty is comparable with the observational 
error bars. Moreover, comparing results in figO shows that 
determining the optical depth hardly depends on the as- 
sumed rotation curve, suggesting indirectly that the result 
is not very sensitive to changes in parameter in the range 
from 200km/s to 239km/s. For comparison, in addition to 
the rotation curves A and B, we consider the third rota- 
tion curve, which is based on the terminal rotation curve 
data from (|McMilladl2LTllh with the assumption of a higher 
value Oo = 239km/s. 



2 In the simulation that also takes the bar into account, the full 
range I = ±5° is considered (see section [2. 711 
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Figure 5. Optical depths in disk model with the disk scale height 
h = 325 pc. Points represent the observational data averaged over 
longitude within the range out to £ = ±5°. Top panel: The result 
for the rotation curve A [thick line] and after subtracting the gas 
contribution [thin line]. Bottom panel: Comparison of optical 
depths for different model rotation curves: A [solid line] and B 
[dashed line] and, in addition, the rotation curve based on the 
McMillan data [dotted line] (without gas subtraction and with 
Ro = 8.29kpc preferred by McMillan). 



2.4 The influence of the solar Galactocentric 
distance R a 

In our calculations we adopted the frequently used value of 
the s olar Galac t ocent ric distance R — 8.0 kpc, the same 
as m lAvedisoval l|2005h . where a review of other estimations 
can be found. In modeling microlensing in the Galaxy, other 
values are sometimes preferred. In the Galactic model pre- 
sented by the EROS collaboration (De rue fc at al|[l999l ). the 
value recommended by the IAU R a = 8.5 kpc is used. Be- 
cause the integrand in the formula defining the optical depth 
((3| is small near the Sun's position, the choice of R can- 
not alter the results significantly, however. When R a is in- 
creased from Ro = 8.0 kpc to R a = 8.5 kpc in the integral 
Q), there is a corresponding change in the optical depth of 
At = +0.1 x 10~ 6 at b = 1° and only At = -0.04 x 10~ 6 
at b = 5°. One can also find a lower value for R a of 
7.62 kpc (|Eisenhauer et all 120051 ). Replacing R a = 8.0 kpc 
with Ro — 7.62 kpc in the integral ((3| gives rise to a change 
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Figure 6. A Monte Carlo simulation of the optical depth (and 
standard deviation) at various latitudes b [solid circles+bars], 
compared with another model (discussed earlier) assuming the 
sources of light to be located on the symmetry axis [solid line]. 



in the optical depth of Ar = —0.1 x 10~ 6 at b = 1° and only 
At = +0.02 x 10~ 6 at b = 5°. This shows that the optical 
depth is not sensitive to the changes in the estimate of R a 
within the limits most frequently met in the literature. 

2.5 The influence of the bulge structure 

The disk model predicts a substitute surface density o~(R) 
in the disk plane that accounts for the rotation curve. It 
can be regarded as the column mass density. If the true 
mass distribution is flattened then the description in terms 
of the column mass density should be sufficient for a dy- 
namical model, but an appropriate volume mass reconstruc- 
tion is needed when studying the microlensing. In the pre- 
vious section we assumed that the vertical profile of the 
total volume mass density p(R, z) changed exponentially: 
p(R,z) = p(i?,0)e- |z|/ ' 1 = a(R)e^^ /h /2h. In particular, 
the central bulge can be described in terms of its own column 
mass density, which contributes to the total mass density. 
But the observations of the bulge (for example. iBinnev et al.l 
Il997l ; iBissantz fc Gerhardll2002l ) indicate that a spherically 
symmetric distribution of matter or even a three-axial el- 
lipsoid would better approximate the central part of the 
Galaxy, and close to the center there can be some correction 
to the density profile we use. The stars residing in the cen- 
tral bulge contribute significantly to the optical depth, thus 
the particular geometry of mass distribution may have some 
effect on the accuracy of our result. Therefore, some estima- 
tion of the uncertainty is indispensable in our predictions 
concerning the optical depth. 

To this end, we replace the central part extending out 
to lkpc by a spherically symmetric bulge and check how 
this change in the mass distribution alters the optical depth. 
Since the dynamical mass inside lkpc is comparable in both 
models (the disk model gives Mi k p C = 1.061 x 1O 1O M , 
whereas for the spherical bulge Mitpc = 1.137 x 1O 1O M0), 
we do not expect large discrepancies in their predictions. A 
possible change by a factor of order of unity may result from 
a change in the iso-density surfaces and from a small change 
in mass distribution outside lkpc compensating for any ex- 



cess from the rotation that including a spherically symmetric 
part may have caused in the gravitational potential. 

With the bulge included, the integral for the optical 
depth splits as follows 



s(l - s)p diak (s)ds + / s(l - s)p bu i ge (s)ds 



c 2 r h (b) 
AttGR 2 



The first term is the contribution from the thick disk, and 
the second term is the contribution from the central bulge, 
where pbuige and pdisk are understood as functions of R(s) 
and z(s) along the lines of sight. The integration limits are 
determined by the bulge size Rb = (1 — Sb)Ro = lkpc. The 
density pbuige at the spherical radius r(s) — V R 2 + z 2 — 
R ^(l-s) 2 + s 2 b 2 can be found from the mass function 
M(r) of the bulge 



pbulge(s) = 



dM 



4nr 2 (s) dr 



r = r(s) 



and this M(r) in turn contributes to the rotation curve. 
We assume that the bulge accounts for the measured rota- 
tion curve for radii R < Rb- In this case we have M(r) — 
rv 2 (r)/G, while for R > Rb we assume the contribution 
from the bulge to be Keplerian: Vb(R) = v(Rb)y/Rb/ R- 
As previously, for small &'s the disk volume density is 
Pdisk = (Ti(i?o(l — s))e~ s ^ R °/ h /2h, where cri accounts for 
the remaining part of rotation v 2 (R) — v 2 (R) — v 2 (R) not 
explained by the bulge (a similar de composition of rota - 
tion was applied to another galaxy in lJalocha et al.ll2010l ). 
The decomposition of the rotation curve is plotted in fig[7] 
The required surface density o~i(R) corresponding to vi(R) 
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Figure 7. The decomposition of the measured rotation curve into 
a bulge and a disk contribution. 



ca n be found by app lying the iteration method described 
in lJalocha et al. (|2008l ). The resulting optical depth is plot- 
ted in figlU As expected, the correction from the presence 
of the bulge is so small (of the order of Ar = 0.2 x 10 -6 ) 
that it cannot change the main conclusion of the previous 
paragraph. 
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Figure 8. The optical depth in a model with a spherically sym- 
metric bulge [thin dashed line] and the disk model prediction 
[solid line]. 



2.6 Influence of the vertical structure 

The disk model predicts only the column mass density u(R), 
and its reprojection to a volume mass density p(R, z) de- 
pends on the additional assumption about the adopted den- 
sity fall-off in the vertical direction. It is therefore indis- 
pensable to estimate the influence of this assumption on 
the predicted optical depth. As the starting point, we adopt 
a single exponential vertical fall-off p(R, z) = p(R, 0) e~^ z ^ h 
with the scale height parameter h = 325 pc. The dependence 
on the scale height parameter has been already checked and 
presented in figure [4] It shows that when the effective change 
of the scale height is not large, the results cannot be changed 
significantly, especially for the higher latitudes. In this sec- 
tion, we compare our result with the ones for two other 
vertical distributions of matter. First, we consider a double 
exponential vertical profile: 



Pdisk(R,z) = p{R,0) (l-/3)e 



\z\/h ± 



+ /3e 



The parameters hi — 320 pc for the thinner disk, /12 = 
643 pc for t he thicker di s k, an d the fraction /3 = 0.216 are 
taken from [Gould et all (1 19971 ). L ikewise, we also conside r 
the hyperbolic cosine profile used in lvan Altena et"aL 

1 

pdisk{R,z) = p(R,0) 



cosh(z/h) 



As we show in section I2l9l this profile approximates the ob- 
servations better in the solar neighborhood and improves 
the mass-to-light ratio. The comparison of the model curves 
r(b) for the disk vertical density distribution profiles men- 
tioned above is presented in figure [9] The calculations are 
made after subtracting the gas contribution and by consid- 
ering the spherical bulge component. As we expected, the 
choice of the particular vertical density fall-off could only 
slightly change the results. Replacement of a single expo- 
nential profile by a double exponential profile reduces the 
optical depth approximately by At = 0.2 x 10 -6 . 

Finally, we discuss a more complicated situation. The 
case of a three-axial ellipsoid is beyond the scope of the 
simple axi-symmetric dynamical model of the Galaxy, nev- 
ertheless, a correction to the optical depth connected with 




* OGLE all stars 
A MOAall stars 
□ MACHO all stars 
x OGLE bright stars 
MACHO (1 year) 
■ MACHO bright stars 

EROS {3 years) 

• EROS bright stars 
disk vertical profiles: 

exponential fall-off 

double exponential 

fall-off 
1 profile 



1.5 2 2.5 3 3.5 4 4.5 5 5.5 

b [deg] 

Figure 9. The optical depth calculated after subtracting the gas 
contribution, with the inclusion of the spherical bulge, shown for 
all of the vertical density distributions considered in the text. The 
thin dotted line represents the reference model - a disk with the 
exponential vertical profile with the scale height h = 325 pc. The 
thick dotted line corresponds to the disk with a double exponen- 
tial vertical profile, and the solid line represents the hyperbolic 
cosine vertical profile. 



the presence of the bar should be within an order of magni- 
tude as in the case of the spherically symmetric bulge con- 
sidered earlier. We have seen above that our result is very 
stable against small changes in the mass distribution. Sim- 
ilarly, considering a bar-like structure would only introduce 
higher multipoles to the mass distribution that should not 
affect qualitatively our main conclusion. In the next section 
we study this issue in a more detail. 



2.7 Model uncertainties and the bar issue 

It is plausible that very particular assumptions about the 
dynamical model of the Galaxy could affect the resulting 
optical depth. If the model is based on the rotation curve, 
in preparation of which the axial symmetry is assumed, it 
is natural to assume axisymmetry for the distribution of 
the sources of light. This symmetry is consistent with the 
disk symmetry and with the spherical symmetry, whereas 
it is not consistent with the three-axial symmetry of a bar- 
like structure. Because of this, an axisymmetric model is 
not adequate for estimating the optical depth in a single, 
particular direction. That is why we are using a t(6) ob- 
servable, which has the meaning of the optical depth over 
slices of fixed longitude with latitudes around a particular 
b. Taking the statistical character of r(b) into account, we 
performed a Monte Carlo simulation that mimics the ob- 
servational situation of detecting microlensing events from 
many sources with various longitudes and various distances 
to the observer. Although the longitude of each field is well 
determined, the distance to the particular source of light 
(described by the x parameter) is practically unknown. We 
have chosen a sample of pairs (I, \) f° r which we used a prob- 
ability distribution weighted by the surface mass density a. 
The result for the r(6) calculated for this sample is depicted 
in figure 6. It shows that on average the approximation for 
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which the sources are located on the symmetry axis works 
quite well. 

One of the advantages of the Monte Carlo approach is 
that it enables one to estimate a possible change or uncer- 
tainty in the determination of r that could be introduced 
by deformations of the central part, in particular those due 
to the presence of the bar. As previously, we do not expect 
changes in the main conclusion. This remark is suggested 
by the fact that the dependence of the optical depth on the 
longitude t(1) (within latitudinal slices of longitudes from 
around a particular I) is much we aker than the dependenc e 
on the latitude r(b). According to Hamadac he et all |2006), 
the optical depth changes by about Ar = 1.3 x 1CP 6 in the 
range I £ (—5°, 5°) while the optical depth change in the 
considered range of latitudes is about Ar = 2.7 x 10 -6 . Our 
disk model predicts the change in the optical depth in the 
range I g (-5°, 5°) of about Ar = 1.8 x 10" 6 . Then, the 
correction connected with the inclusion of the bar should be 
less than Ar = 0.5 x 10~ 6 . To estimate the influence of the 
bar on our results, we e mploy the three-axial G2 model de- 
veloped in iDwek" et all l | 1993 1. It gives a volume emissivity 
of the sources of light in the Galactic center: 



MonteCarlo: three-axial bulge model 
MonteCarlo: axial symmetry 
reference model 




Figure 10. A Monte Carlo simulation of the optical depth (and 
the standard deviation) in the axisymmetric case [solid circles 
and bars] and in the presence of a three-axial bar [open circles 
and bars] . The solid line represents the model curve based on the 
integral J2), in which case the sources of light lie on the symmetry 
axis. 



ri(x',y',z') = ??oexp 



+ 



•(4) 



In section [2721 we used the following notation for the carte- 
sian coordinates of the source X® = [x, y, z] — R a [1 — (1 + 
X) cos 6 cos Z, — (1 + x) cosbsinZ , ( 1 + y) sin b], which differs 
from the naming convention used in lDwek et al.l |l995) . If we 
denote by x' , y , z as the cartesian coordinates of their co- 
ordinate system, then the appropriate transformation reads 
as x' — > y, y' — > —x, and z' — > z. 

We applied the formula @ to perform our Monte Carlo 
simulation in which the sources of light are chosen randomly 
with the probability distribution proportional to the volume 
emissivity. This imitates the situation in which the sources 
are distributed within the three-axial bar. We take the pa- 
rameters xq — 2.01 kpc, yo = 0.62 kpc, and zq = 0.44 kpc 
fitted in the 3, 5 /mi band, although the values obtained in 
other bands give roughly the same result. In figure fTUl below. 
we present the comparison of the Monte Carlo simulations 
assuming the axial symmetry and one that considers the G2 
three-axial model. For this three-axial model, an extended 
range for longitudes is assumed j I j< 0.085 (that is, ±5°), 
which covers all t he fields used in the pre paration of data 
(see figure 14 in lHamadache et al.l [2006) . In the central 
part of the bulge (say, ±1° in the longitude), there is practi- 
cally no difference between the predictions of the three-axial 
model and the spherically symmetric bulge models. One 
can see that the results for the sources distributed within 
the three-axial bulge are also consistent with the main disk 
model prediction. The correction due to the bar, combined 
with the dependence on the choice of the field range, gives 
no more than Ar = 0.4 x 10 -6 . We neglect that the non- 
axisymmetric component affects the dynamics of the inner 
part of the Galaxy a bit. However, this correction also cannot 
be large, because the photometrically det ermined bulge mass 
in the G2 model Af G2 = 1.3 x 10 10 M Q IDwek et all (|l995l ) 
is comparable to our estimation Mik pc = 1.137 x 10 10 Mq. 



2.8 The hidden mass estimation 

To estimate the amount of the hidden mass we must take 
the r(b) prediction in the thin disk model framework into 
account after subtraction of the gas contribution, using rota- 
tion curve B (which gives a more massive disk than rotation 
curve A), considering the spherical bulge component and the 
disk double-exponential vertical fall-off. The result can be 
parameterized as follows: r(b) = ciexp(— C2 [ b |) X 10 -6 with 
coefficients ci = 4.23±0.65 and c 2 = 0.228±0.056. Not only 
the neutral and molecular hydrogen and helium, but also 
other components of the total dynamical mass inside the so- 
lar orbit of Mtot = 6.17x 10 10 M Q contribute to this result. A 
fraction of Mtot could consist of matter that is neither in the 
form of compact objects nor in the form of hydrogen and he- 
lium clouds, which have already been subtracted. Because it 
would not contribute to the measured optical depth, we call 
it the hidden mass. It overestimates the model optical depth 
r(b) by a factor Thidden{b) , which can be obtained by insert- 
ing the volume mass density of the hidden mass phidden(R, z) 
into the integral @. Unfortunately, the Phidden(R, z) profile 
is not known. Since this profile appears in the integrand of 
the relevant formula, there are various possible profiles of 
Phidden(R, z) that could lead to the optical depth correction 
Th t dden(b) giving T(b)-T hid de n (b) as the best fit to the obser- 
vational optical-depth data points. For the simplicity of the 
estimation, we may assume that the hidden mass density is 
a constant fraction of the dynamical mass density (without 
gas) phidden(R, z) = f p(R, z) . Then, the correction factor 
is Thtddenib) = f t (b) , and the observed optical depth is ex- 
pected to be (1 — /)r(b). We find the fraction / by fitting 
the formula (1 — /) ciexp(— c 2 b |) x 10 -6 to the bright 
stars r observations, by means of the least squares method 
weighted by the observational uncertainties. The results of 
the fitting procedure are depicted in figure 1111 Taking the 
OGLE, MACHO, and EROS data into account, we obtain 
/ = —0.139 ± 0.059. Because this fraction is negative, that 
would mean that an additional dynamical mass of about 
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OGLE bright stars 
MACHO bright stars 
EROS bright stars 
tit with OGLE data 
fit without OGLE data 




b |dcg| 



Figure 11. The (1 - /) ciexp(-c 2 | 6 |) X 1CT 6 fit to the 
OGLE, MACHO, and EROS bright stars optical depth data 
[solid line] and a similar fit to the narrowed data of MACHO 
and EROS bright stars data [dashed line]. The filled area corre- 
sponds to the standard deviation range of the Monte Carlo sim- 
ulation (including bar) described in the previous section. The 
ciexp(— C2 | b |) X 10 -6 curve is not plotted for clarity. 
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Table 2. The expectation value, the variance, and the test statis- 
tics z;, in the case where all data are considered 



I f M tot |= (0.86 ± 0.36) x 1O 1O M is needed. The data 
points collected by the OGLE collaboration lie above the 
points collected by other groups and have a bit larger un- 
certainties. If we restricted the analysis to the MACHO and 
EROS data, we would get / = 0.018±0.067, which translates 
to | f Mtot |= (0.11±0.41)xl0 10 M Q of the hidden mass. One 
should regard all this as rough estimates. The two fits that 
we present are only used to find the expectation value of the 
optical depth E[r(b)] by means of the least squares method. 
On the other hand, we know what the expected variance of 
the optical depth predicted by the thin disk model is, be- 
cause we have performed the Monte Carlo simulation, which 
gives us the uncertainties ATMonteCario{b) ■ They could differ 
from the fit residuals. The variance assigned to our model 
should be var[r(6)] = (ArMonteCario(b)) 2 ■ Instead of check- 
ing the goodness of the fit, one can compare the model 
distribution {E[r(b)], var[r(fo)]} with the observational one 
{f(b), a 2 (b)}. The average f(b) can be taken within some 
narrow band in the latitude 6±0.2deg, and the a(b) is the ob- 
servational uncertainty within the band. To check the statis- 
tical hypothesis that these two distributions have the same 
expectation values E[r(b)] = f(b), one should compute the 
quantity: 

^ = E[r(b)]-f(b) 

b / var[T(6)] . <t2(6) ' 

V n n b 

Here n = 10000 denotes the number of points generated in 
the Monte Carlo simulation for the fixed latitude b, and rib is 
the number of the observational points within the band 6 ± 
0.2deg. The variable Zb has the normal distribution Af(0, 1). 
Below, in Tables Q] and [21 the appropriate values are listed 
in the case without the OGLE data and in the case in which 
all the data are considered. 

In both cases the test statistics Zb have values relatively 
close to zero. On a reasonable confidence level of p = 0.01 
one cannot reject the null hypothesis that the model distri- 



bution {J?[t(&)], var[r(6)]} and the observational distribu- 
tion {f(b),a 2 (b)} have the same expectation values. By tak- 
ing the highest value Zj = 1.3, one could reject this hypoth- 
esis on the confidence level p = 0.19; however, this would 
lead to a very high probability p — 0.19 of the statistical 
error of the first kind (of the rejection of the true null hy- 
pothesis). With a smaller number n, the agreement is even 
better (because for n — 10000 the term var[r(6)]/n can be 
neglected) . 

In other words, almost all observational data lie close to 
the expectation value E[r(b)], within the ±AMonteCario(b) 
range, which is depicted in figure ITTI Since the fits we have 
obtained are close to the value / = 0, and the observational 
data are consistent with the Monte Carlo simulation, we 
are led to the conclusion that all of the dynamical mass 
corresponding to the rotation curve in the disk model (and 
after subtraction of gas) is consistent with the optical depth. 
The amount of the hidden mass or an additional dynamical 
mass that could be needed to account for the optical depth 
observation, is negligible in the < R < R a region. 

2.9 The mass-to-light ratio in the solar 
neighborhood in the disk model 

As a cross-check for the dynamical model of the Galaxy one 
can ask about its prediction of the light-to-mass ratio in the 
solar neighborhood. To determine the profile of the local 
mass-to-light ratio one needs a brightness profile. It is dif- 
ficult to determine the brightness profile for Galaxy in the 
form of a single curve. There were only attempts to deter- 
mine the overall brightness in the solar vicinity by direct 
counting of nearby stars. Such brightness estimate in our 
neighborhood is quite accurate (but questionable as glob- 
ally repres entative for t he Ga laxy at the same radial dis- 
tance). In iFlvnn et alj (|2006l ) the local I-band brightness 
is 29.54Lq/pc 2 . It is determined based on star counts out 
to 200pc above the Galactic midplane. In the same volume 
there is 13.2M /pc 2 of hydrogen (HI + HII + H 2 ). Given 
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these data, we obtain IIOMq /pc 2 or 14QMq /pc 2 in the disk 
model for the column mass density corresponding to rotation 
curves A or B, respectively. Subtracting four thirds of hydro- 
gen (thus including He from the bariogenesis) , because the 
gas is nonluminous in this band, we obtain M/Li — 3.13 or 
M/Li = 4.14, respectively, in the solar vicinity. Although 
overestimated, these values are low (and could be still re- 
duced by including more distant objects, above z = 200pc, 
which have not been included in the counting). 

For a better estimate of the mass-to-light ratio in the so- 
lar neighborhood we can test other models of dis t ribut ion of 
stars above the midplane. In Ivan Altena et all ((2004), the 
profile ansatz 1/ cosh(z//i) was assumed as approximating 
the observations better than the usual exponential falloff. 
Assuming our total (column) surface density of IIOMq/pc 2 
and h = 330pc we obtain p = O.IIMq/pc 3 in the solar 
vicinity, so that the integrated density out to h — 200pc 
gives 42Mq/pc 2 and the corresponding M/Li = 0.83 or 
M/Li = 1.1 for rotation curves A or B, respectively. With 
exponential vertical falloff e~' z '/ h (with h = 350pc), we 
would have M/Li = 0.95 or M/Li = 1.47 for rotation 
curves A or B, respectively. These are quite reasonable and 
low local mass-to-light ratios, which seem to be compatible 
with the recent findings from the star k inematics in the solar 
neighborhood l|Moni Bidin et alj|20l3 ). 



3 SUMMARY AND CONCLUSIONS 

The microlensing method determines the optical depth ob- 
servable. The observable can be modeled provided the vol- 
ume mass density of compact objects is known. We assumed 
that the total volume mass density has an exponential ver- 
tical falloff and that the associated column mass density 
accounts for the tangential component of rotation in the 
thin disk model approximation. The structure of the central 
bulge and the presence of the bulge do not alter the re- 
sults significantly, so our results are also insensitive to slight 
changes in the rotation curves such as the choice of the pa- 
rameter O.) Having subtracted the gas contribution from the 
volume mass density obtained in this way, we modeled the 
optical depth and then compared the model with the optical 
depth measurements. It follows from the analysis that the 
amount of mass seen through gravitational microlensing in 
the region < R < R a is consistent with our mass model. So 
that our model of the optical depth conform more to realis- 
tic measurements, we carried out a Monte Carlo simulation. 
It proved to mimic the measurements quite satisfactorily, 
since on the reasonable confidence level one could not re- 
ject the hypothesis that its expectation value equaled the 
expectation value of the observational data. We estimated 
the limits for the amount of nonbaryonic dark matter in our 
model inside radius R to be (0.11±0.4) x 10 10 M Q . Since mi- 
crolensing only detects compact objects, the result suggests 
that the nonbaryonic mass component may be negligible in 
this region because the mass of gas and that of compact ob- 
jects seen through microlensing suffice to account for Galaxy 
rotation in this region. 
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APPENDIX A: THE INTEGRAL 
TRANSFORMS APPEARING IN THE DISK 
MODEL 

The relation {TJ used in the text can b e obtained from an 
equivalent integral expression derived in (jBratek et apj ool) 



2 [ K( X /R) _R E( X /R) ] 



To this end one can use the identities 2E (k) = 
(1 + k) E ( Mfe ) + (1 - k) K ifi k ) and (1 + k) K (k) = K (p k ) 
with fit = §+f> applied separately for \ < R (with fc = 
and f° r R < X (with fc = R/x)- The inverse rela- 
tion mapping the surface density to the rotation law can be 
found directly from the definition of the gravitational po- 
tential, without the usual intermediate step through Henkel 
transform^] (note the sign difference in {TJ) 



<*> = 2G 



R 



1 + x 1 



a(Rx) xdx. (Al) 



The two integrals {TJ and (|A1[) are quite symmetric, the 
inversion transform x — > l/x leaving fi x invariant links the 
integral kernels of the two integrals. We see that two physical 
quantities, the disk surface density cr(R) and the centrifugal 

acceleration on a circular orbit " x 1 ' in the disk plane, are 
integral transforms of each other. It would be interesting to 
investigate this striking property of the disk model in more 
detail in the context of the theory of integral transforms. 



3 The gravitational potential due to axisymmetric thin disk with 
surface mass density a in cylindrical coordinates is 

oo 

i (Rs)= i G f . *" (x)dx .. r d ° 

j ^{r+xy+z 1 Jo ,/r 



> ' \ 1 - , | 4 ^i '1 sin2 ( a ) 



and the definition of elliptic function K can be used. Differenti- 
ation with respect to R and taking the limit z — > (the velocity 
on circular orbits in the disk plane is v 2 (R) = Rdn^ (R, 0)) gives 
the desired result. 



